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I. INTRODUCTION 


Military airborne tactical weapon systems have been 
designed and implemented primarily to aid the aircrew in 
performing their mission with accuracy and speed. This 
usually means that a shorter time is spent over target 
which increases the survivability of both the weapon 
systems platform and the aircrew. 

The system presently employed on board the A6E, one 
of the Navy's attack aircraft, utilizes an IBM 4 PI series 
mini-computer to perform two major functions. 

1. Navigation. 

2. Solution of the ballistics problem. 

In addition, several other related functions are performed 
by the system. 

3. It provides steering commands. 

4. It provides real-time display of sensor information. 

5. It provides release pulses to the weapon at the 


appropriate time. 


A. FEASIBILITY STUDY 


This thesis will attempt to prove the feasibility of 


implementing an airborne tactical weapons system using 
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micro-computers. Two important questions must be answered 
in order to establish the feasibility of using a micro- 
computer in a tactical weapons system: 

1. Is the micro-computer accurate enough? 

2. Is the micro-computer fast enough? 
The accuracy problem was approached by executing the ballis- 
tics algorithm on an IBM system 360 using a 21 to 24 bit 
mantissa in the standard floating point format and comparing 
the results of the same algorithm executed on an INTEL 8080 
micro-processor with a 16 bit floating point mantissa. The 
question of speed was answered by aeeeueane the ballistics 
algorithm for numerous weapon types and initial conditions 


and observing the elapsed clock time. 


B. VERIFICATION 

The second aspect of this thesis was to verify that the 
ballistics ee eee to published tables as 
well as experimental data. The NAVAIR 01-1C-1T-1 ballistics 
tables were used to compare time of fall and down range 
travel against the results of the FORTRAN version of the 
ballistics algorithm. A total of 1813 different initial 
release conditions were examined spanning 18 different 
weapon types, various dive angles from +10 degrees to -60 


degrees, altitudes from 500 feet to 15,000 feet, and air- 


speeds from 300 knots to 650 knots. 
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The source of experimental data was a set of data cards 
recorded by the bombardiers of an A-6E Squadron over a 
period of one year. The data for each bomb drop consists 
of 24 various computer readouts at the instant the weapon 
is released from the aircraft as well as the hit coordinates 
of the weapon. This information is used to determine the 
initial conditions for the ballistics equations. Unfortun- 
ately a critical parameter, the dive angle, was recorded 
only to the nearest degree. Consequently an error analysis 
was conducted to determine the maximum error which could be 


expected from the rounding of the output data. 


C. PRESENTATION OF THE THESIS 

Chapter II explains the organization of a distributed 
micro-computer airborne tactical weapons system and dis- 
cusses how the output of each subsystem is integrated with 
the entire system. The ballistics problem and the deriva- 
tion of the differential equations used to-describe the 
mathematical model are discussed in detail in Chapter III. 
Since no analytical solution exists for these equations, a 
simplified version of the model is solved analytically to 
gain insight. 

Chapter IV explores the feasibility of using micro- 


processors in an airborne tactical weapons system. 
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Aspects of accuracy and speed are examined. Factors in 
attempting verification of the ballistics algorithm and 
experimental a aren the NAVAIR 01-1C-1T-1 ballistics 
tables are re in Chapter V. 

Chapter VI presents the results of this thesis. The 
conclusions and recommendations of the author concerning this 
thesis are written in Chapter VII. 

Appendix A contains the output from the FORTRAN program 
comparing the FORTRAN and PLM versions of the ballistics 
algorithm. Appendix B presents the same comparison between 
the FORTRAN version of the ballistics algorithm and the 
ballistics tables. Appendix C is a listing of the experi- 
mental data as they were recorded from the cockpit readouts 
of various A-6E aircraft. Appendices D and E are the results 
of the experimental delivery data compared against the 
FORTRAN version of the ballistics algorithm (approximating 
the ballistics tables), using two sets of drag coefficients. 

The ballistics algorithm used in this thesis is the 
SIGMA version of the BOEING ALGORITHM modified by the Naval 
Weapons Laboratory at China Lake, California. The algorithm 
was further modified at the Naval Postgraduate School in 
Monterey, California for eventual implementation on the 
ballistics processor of the multiple micro-processor tactical 


weapons system. 
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II. BACKGROUND 


To provide a better understanding of the role the 
ballistics processor plays in the multiple processor 
system, a brief overview of the proposed system will be 
discussed. The computer system is composed of three micro- 
processors: a navigation computer, a ballistics computer, 
and an executive computer. Each machine is dedicated to 
the process to which it is assigned instead of sharing 
resources of a single processor as in the present opera- 


tional systems. 


A. THE NAVIGATION COMPUTER 

The navigation computer is a basic element to all 
tactical systems. In present operational systems, the 
navigation program is executed periodically to update the 
present position by the change in position since the last 
time increment. The navigation computer utilizes input 
from four major sensor instruments as its primary source 
of information. 

l. The Inertial Navigation System. 

2. The Doppler Radar. 

3. The Air Data Computer. 


4. The Search Radar. 
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The inertial navigation system provides heading 
(azimuth), attitude (roll and elevation), and velocity 
increments in the X, Y, and Z directions. The doppler 
radar is a velocity sensor that utilizes the doppler shift 
principle = measure ground track speed and drift angle. 
The ground speed and drift angle derived from the doppler 
radar along with the true heading from the inertial naviga- 
tion system are used to calculate the direction and magni- 
tude of the wind. The air data computer uses the ambient 
static pressure and ram air pressure from the pitot tube 
to calculate WeFeected static pressure, pressure altitude, 
and mach number. Outputs from the air data computer are 
used to damp the raw velocities from the inertial system. 
The search radar provides target aximuth, range, and eleva- 
tion signals to the navigation computer. The search radar 
elevation along with the aircraft elevation from the inertial 
system are used to measure the radar depression angle (look 
down angle to target from flight path vector). The depres- 
Sion angle and search radar slant range are used to compute 
ground distance to the target independent of altitude. 

As the navigation computer calculates new incremental 
distances for each time increment, the present position 
is continuously updated. After each update, the executive 


computer is interrupted and the current value is passed to 


ies: 





it. If the navigation computer is functional, the aircrew 
will have current present position information independent 


of the status of the executive computer. 


B. THE BALLISTICS COMPUTER 

The ballistics computer is provided with the most current 
estimates of position and velocity of the aircraft. It is 
also provided with the weapon type selected by the aircrew. 
The ballistics algorithm computes the down range travel 
and time of fall for the weapon based on the airspeed, dive 
angle and altitude of the aircraft. The ballistics algorithm 


will be discussed in depth in Chapter III. 


C. THE EXECUTIVE COMPUTER 

The executive computer displays the information generated 
by both the ballistics and navigation computers. The execu- 
tive computer is also responsible for issuing steering 
commands to the autopilot and firing pulses to the weapons 
release mechanism. The most important task the executive 
Bibate: performs is to extrapolate a predicted weapon release 
point based on a time history of the position and velocity 
of the aircraft. Thus, even though the other two computers 


can operate independently, the aircrew cannot make a computer 


delivery without the executive computer. 
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D. THE COMMUNICATIONS SCHEME 

Because a5 the inherent hierarchy between the computers 
involved, the master-slave type of multiprocessing is the 
most suitable and simple form to implement. The navigation 
and the ballistics computers act as peripheral devices of 
the executive computer, resulting in a one way interrupting 
scheme. The only computer which has to be interrupted is 
the executive computer. The navigation and ballistics 
computers are the dedicated slaves which asynchronously 


interrupt the executive computer. 


Ly 











INERTIAL 


NAVIGATION 
SYSTEM 





BALLISTICS COMPUTER 


FIGURE 1. ORGANIZATION OF MULTIPLE PROCESSOR AIRBORNE 
TACTICAL WEAPONS SYSTEM 
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III. BALLISTICS PROBLEM 


Since the earliest days of aerial warfare, She, recess of 
the ballistics problem has been to drop a bomb from an air- 
borne weapon platform and to consistently hit a target. The 
problem still exists today; however, with the aid of a 
computer the problem can be quickly and accurately solved. 
This necessitates the development of a mathematical model 
which approximates the actual path a weapon travels through 
Space. Statistically, projectiles have been shown to follow 
predictable paths which behave very much like freely falling 


bodies described by Isaac Newton. 


A. THE MATHEMATICAL MODEL 

The ballistics problem can be described as a body 
falling through space according to Newton's second law of 
motion. The mathematical model or equation governing the 
trajectory of a ballistics projectile is a second order 


differential equation: 


ma = mg - c [viv (1) 


where: m = mass of the body 


acceleration vector 


4) 
ll 


acceleration due to gravity 


the) 
Il 
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e drag coefficient 


v = velocity vector 
The term e( Iv|v) is the drag due to the air resistance of 
the body and is proportional to the square of the velocity. 
At the time the projectile is released from the aircraft 
it will have an initial position and velocity. The differ- 
ential equation together with the initial conditions 
uniquely determines the trajectory of the projectile. 
This particular mathematical model was chosen for two 
reasons: 
l. It approximates reality accurately. 
2. The Navy publishes range and time of fall 
information in tables (NAVAIR 01-1C-1T-1) 
which uses this model. 
The latter is the primary reason for using the ballistics 
tables as a standard or guideline in the verification of 


any new mathematical models or program implementations. 


BRB. BASIC ASSUMPTIONS 

Unfortunately there is no known analytical solution for 
this set of differential equations. However, a simpler 
problem does have an analytical solution. Thus, for the 
sake of simplicity and to aid in the discussion of this 


solution, let us first assume: 
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1. Level, non-accelerated flight. 
2. Bombs are not ejected from the bomb rack. 
3. No forward firing ordnance. 
4. Drag/wind resistance is negligible. 
5. Time of release is time = 0. 
In addition, this discussion also assumed: 
6. The Earth is flat and non-rotating. 


7. The gravitational attraction, g, is constant. 


C. DERIVATION OF DOWN RANGE TRAVEL 
Rewriting equation (1) with respect to the time deriva- 


tives of the position vector, say u, the equation becomes: 


mu = mg . Cc jul u (2) 


By letting u = (x,y), equation (2) is transformed into a 
system of two equations describing motion in a two dimen- 


sional coordinate system. 

mx |. Pa | x (3) 
‘eo 22 72 : 

my = -mg + Vx ty | (4) 


Initial conditions for the system are derived from the re- 


lease position and velocity vectors of the aircraft: 


72), 





Position Velocity 
x(tp) = r, (to) x(to) = r,(to) 
y(ty) = 1. (to) y(to) = (ty) 


As a result 


negligible and the aircraft's velocity and 


mine the initial conditions, equations (3) 


simplified to 


mx = 0 
my = -mg 


The solution of 
by the constant 


and solving for 


of the assumption that air resistance is 
position deter- 


and (4) are 


initial conditions 


x(0) = r,(0) x(0) = 15 ACY, 

KO) = PG 0) = r(0 

y(0) r,( ) y(0) = r,(0) 
these two equations is obtained by dividing 


m, integrating twice with respect to time, 


the constants of integration. 


fsce dt = fo abe = x(t) = Cy 
[ice dt = fe, de > x(t) =et te, 
when t = 0: c, = x(0) = r,,(0) 
c, = x(0) = r_(0) 
therefore, 
x(t) = r(0)t = r_(0) (5) 


where: 
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x(t) = total down range travel after release. 
r (0) = initial displacement x(0) of weapon 
at the time of release, usually 0. 


r.(0)t = down range travel due to initial velocity. 


. 7 
[rs dt > y(t) 


[ve Eat 
7 2 
fect, dt > y(t) = -ggt° + c,t + Co 


[uo dt 


216 2 Cy 


when t = 0: 25 = y(0) = r (0) 
Co = y(0) = r,,(0) 
therefore, 
y(t) = -ket™ + E,(o)t + x (0) (6) 
where: 


y(t) = the height above the ground at any time t. 

r,,(0) = initial altitude y(0) of weapon at time of release. 

ry(0)t = altitude loss/gain due to aircraft's initial 
vertical velocity. 


? 
-3gt = altitude lost due to gravity. 


If the time from weapon release to weapon impact, or 
time of fall, is known, the ballistics problem is reduced 
to determining the down range travel, x(t*), given its 


initial-velocity, r,(0) and the time of fall, t*. Time of 
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fall can be found by setting equation (6) to 0 and solving 


for the positive root of t. 
-ket* + + (0)t + xr (0) = 0 
yi ye 


Using the quadratic formula: 


t* fs - b2 = ie | 2a 


where: 
a = -g/2 
b = ry (0) 
c= el 
ce = Fy(0) + YCr,(0))? + 2¢4,(0) 


g 
By substituting time of fall, t*, into equation (5) 
x(t*) = r (0)t* + r (0) 
x x 


down range travel, DRT, can be calculated. 


= (0) + V/(z.(0))* + 2er_(O 
pat = x(t) = 2,(0) xc * Vityt0)) * Paty). 60) 


& 


Since the coordinate system is arbitrarily placed, 
assume initial displacement, r (0) , to be zero at the time 


of release. Also, since level non-accelerated delivery is 
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assumed, initial vertical velocity, r,(0), is zero. The 
expression for down range travel now becomes: 


DRT = r (0) ary (0) 


a. 
In reality the problem must take into account the con- 
straints placed on the problem earlier. Wind, drag, 
vertical velocity, and non-level delivery prarmeters make 
the solution more difficult. Solution in this case is 
accomplished numerically using a second order Runge-Kutta 


scheme on a digital computer. 


25 





IV. MICRO-COMPUTERS IN AIRBORNE TACTICAL WEAPON SYSTEMS 


The question of whether or not micro-computers are 
feasible in an airborne tactical weapon systems environment 
was approached by first translating the ballistics algorithm 
into two high level languages, FORTRAN -and PIM. (See Ref. 
3). The FORTRAN version was executed on the IBM system 360 
and the PIM version was executed on the INTEL 8080 micro- 
processor. Thus, if the same algorithm is executed on two 
different machines, the INTEL 8080, whose floating point 
mantissa has 16 bits and the IBM 360, whose floating point 
mantissa has 21-24 bits, then the differences in results 
can only be attributed to the difference in the precision 


of the two machines. 


A. THE QUESTION OF ACCURACY 

The solution of the ballistics problem requires solving 
a set of four differential equations numerically on a digi- 
tal machine. This necessitates numerous arithmetic opera- 
tions, including multiplication and division. 

Since micro-computers presently lack hardware multiply 
and divide functions, a software package capable of perform- 
ing floating point arithmetic operations must be used. 


However, this requires a significant amount of additional 
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computing time in the solution of the ballistics problem. 

As a result an alternate design in computer architecture 

was explored which utilizes three INTEL 8080 micro-processors 
instead of one, as in the mini-computer systems. In this 
multiple micro-processor system, each processor is dedicated 
to each of the primary functions of the system: executive, 
navigation, and ballistics computations. 

The floating point mathematical package used in conjunc- 
tion with the ballistics micro-processor uses a three byte 
binary representation with a 16 bit mantissa and an 8 bit 
exponent. The mantissa is left justified so that the most 
significant bit is always on and need not be stored, giving 
a full 16 bits of precision. The exponent is expressed as 
a power of 2 where the most significant bit serves as the 
sign bit. This three byte number is used instead of a con- 
ventional four byte scheme in order to reduce the time needed 
to perform the calculations necessary to solve the ballistics 
problem. 

IBM's system 360 computer utilizes a four byte floating 
point hexadecimal number with a 21 to 24 bit mantissa and 
an 8 bit exponent. Since this method requires the first 
byte to have the value 1 to F (hexadecimal), the precision 
of the mantissa can vary between 21 bits (when the leading 


three bits are zero) and 24 bits. The exponent is expressed 
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as a power of 16 and also contains the sign bit. An 
obvious advantage the IBM 360 has is that it does floating 
point arithmetic in hardware which makes it approximately 
two orders of magnitude faster than the software version. 

The INTEL 8080 micro-computer uses 8 bit operations and 
has the option of using double precision (16 bit) operations. 
The double precision operation permits multiplications and 
division to be performed as sequences of 16 bit additions 
and subtractions. Then, if the mantissa is kept left 
justified (16 bit precision), the double precision feature 
can be used to maintain 16 bits of precision throughout the 
calculation. If more precision is desired, such as 21 or 
24 bits of precision, a quantum jump in execution time can 
be expected because of the additional computer cycles 
required. 

It is the intention of this thesis to show that the 
loss of one byte of precision will not significantly affect 
the results of the ballistics solution. The accuracy with 
which a weapon is delivered depends greatly upon the 
accuracy and precision of the sensor supplied information. 
The loss of one byte of precision only affects the sixth 
most significant digit which is an order of magnitude more 
precise than most of the input sensors on board attack 


aircraft. This then is the motivation to compare the 
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results of a FORTRAN ballistics algorithm executed on the 
IBM system 360 with the same algorithm translated into PLM 


and executed on the INTEL 8080 micro-computer. 


B. THE QUESTION OF SPEED 

The second area of interest is the question of speed. 
At the time this thesis was written, an LSI (large scale 
integration) "chip" existed which could perform the floating 
point multiply and divide operations at the cost of $270. 
However, due to budgetary constraints this equipment was not 
readily available for experimentation. The hardware float- 
ing point "chip'’ can execute approximately 100 times faster 
than the software floating point package. For example, a 
multiply operation in the software package takes approxi- 
mately 600 microseconds to execute whereas the hardware 
package executes a multiplication in 6 microseconds. A 
hardware multiply and divide operator was also developed 
and constructed as a micro-computer course project at the 
Naval Postgraduate School and was demonstrated to function 
at 60 microseconds. 

By interrupting the program during execution and record- 
ing the location of the program address register, it was 
determined that the ballistics program spends about 92% of 


its execution time in the floating point package. According 
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to Jupin (Ref. 3) the execution time of each solution was 
proportional to the calculated time of fall. These results 
were confirmed 2 executing nearly 1800 separate calculations. 
The time to calculate the predicted release point proved to 
be Bbout 10% of the calculated time of fall. Linhares 

(Ref. 2) was able to show that the ballistics algorithm was 
fast enough for certain initial conditions, however for high 
airspeeds and low altitude release conditions his extrapola- 
tion technique was not usable. 

A substantial amount of the time the program spends in 
the floating point package (about 92%) is spent in the 
multiply and divide procedures. Using either the commercially 
produced "chip" or the locally constructed hardware multiply 
and divide operator, a significant reduction in execution 
time would result. Although this thesis will not answer the 
question of speed with an unqualified yes, it supports the 


finding that the ballistics processor is fast enough. 
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TABLE 1. 


NUMBER OF PER CENT OF 
TOTAL INTERRUPTS TIME PROGRAM 
NUMBER OF PROGRAM WAS IN/|SPENT IN 
INTERRUPTS FLOATING POINT} FLOATING POINT 
ATTEMPTED 


PACKAGE PACKAGE 


AMOUNT OF TIME BALLISTICS PROGRAM SPENDS IN 
THE FLOATING POINT PACKAGE 
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V. VERIFICATION OF DOWN RANGE TRAVEL 


Two separate verifications of down range travel were 
made, using the ballistics tables as a "standard." First, 
the FORTRAN and PIM versions of the mathematical model, 
previously discussed, were tested against the ballistics 
tables. Second, observed data was compared against the 
FORTRAN version for accuracy in time of fall and down 
range travel. The FORTRAN version was executed on the 
IBM 360 (32 bit machine), while the PLM version was executed 
on the INTEL 8080 micro-computer utilizing a 24 bit floating 


point mathematical package on an 8 bit machine. 


A. THE BALLISTICS ALGORITHM. - FORTRAN VS. PIM 

A straight forward comparison between the FORTRAN version 
and the PLM version of the ballistics algorithm was made 
contrasting the down range travel and time of fall. An 
input/output interface was written to the PIM program so 
that data could be read from a floppy disk and the results 
written onto the same device. The floating point package 
was also modified in order to execute on the INTEL 8080 
and a logic error in the multiply procedure was corrected. 
The FORTRAN program (Ref. 3) was virtually unchanged, 


however a statement was added to the TRAJ subroutine to 
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patch a logic error affecting the second stage trajectory 


calculation. 


B. OBSERVED DATA VS. BALLISTICS TABLES \ 
Since the ballistics tables are considered a "standard" 
against which various types of ballistics results are com- 
pared for validity, a comparison between observed results, 
obtained from the A6-E experimental data, and the ballistics 
Fables (NAVAIR 01-1C-1T-1) was desired to establish a 
correlation between the two. However, several problems 
were encountered in making the comparison. 
1. Predicted Down Range Travel 

The initial conditions of the observed data are 
not compatible with those of the ballistics tables. The 
observed data has initial conditions composed of various 
dive angles, altitudes, and airspeeds, whereas the ballis- 
tics tables' initial conditions are multiples of 50 and 
100 for altitude and airspeed, and zero for dive angle 
(only considering level delivery). Two apparent solutions 
to this difficulty are (1) the error sensitivity tables 
and (2) interpolation of the ballistics tables. However, 
both have disadvantages. 

The error sensitivity tables failed to help because 


the corrections are based on maintaining a constant sight 
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line (mil setting) rather than keeping down range travel 
constant. For example, in a level delivery situation the 
error sensitivities table assumes a constant sight line 
and varies the down range travel by affecting corrections 
to the altitude. WCE Le che eueueade is higher 
than planned and the sight line is maintained, the hit 
will be short. However, the problem using observed data 
requires a constant down range travel to target, applying 
correction for altitude and airspeed. Therefore, the 
error sensitivity tables could not effectively be used in 


this case. 
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The ballistics tables are not linear. Second and 
even third order interpolation would be required to deter- 
mine down range travel on anything other than a cardinal 
altitude or airspeed. This process is extremely tedious 
and time consuming and was abandoned as impractical. 

As a result, an accurate substitute for the ballis- 
‘tics table value for down range travel was sought. Using 
the ballistics tables for the MK-76 mod-5 practice bomb, 

96 data points were selected from all dive angles, air- 
Breeds, altitudes, and tested against the FORTRAN ballis- 
tics program. Time of fall and down range travel compared 
favorably with a mean error of 0.23% for down range travel 
and 0.11% for time of fall. This was considered accurate 
enough be be used as an approximation for the ballistics 
tables' value for down range travel. 

2. Wind Effect 

The wind plays an important and rather subtle role 
in determining the total down range travel. First, assuming 
a no wind condition, the down range travel, DRT, is computed 
in the direction of the true heading. Since mo wind is 
present to deflect the projectile, the true heading and 
ground track will coincide and the down range travel along 


the true heading and ground track will be equal. 
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Now consider the wind. If a coordinate system 
(x,y,z) is introduced such that the air does not move with 
respect to the coordinate system, then an aircraft in this 
air mass has velocity with respect to the air mass and its 
heading will be the true heading. Because the air mass 
moves with respect to the ground with the wind velocity, 
w = (w,, Wy » w,), the Earth fixed coordinate system 


(x',y',z') is related to(,y,z) by the following equations: 


| = 

x =x Fo wy,t 
'sey twit 

y y y 
t = 

Z Z twit 


Considering only the horizontal wind, the vertical 
wind, wz, becomes zero. Down range travel can now be deter- 
mined given an initial altitude, z, and air velocity, V, 
and will be in the direction of the true heading. To 
determine the point (x',y') in the Earth fixed coordinate 


System at which the weapon lands requires: 


*% 
ll 


(GhO Lf) mie eee (tof) 


“< 
Ml 


y(tof) ae * (tof) 


For example, given the same initial conditions, down range 
travel will be computed, as before, in the direction of the 


true heading. However the aircraft will drift with the 


37 





wind and will actually move across the ground on a differ- 
ent heading, called ground track. Since DRT is computed. 
along true heading, the projection of DRT onto the ground 
track is that distance the projectile will travel due to 


the initial conditions alone. 
PROJECTED DRT = cos(TH-GT) * DRT 


The down range component (along the ground track) of DRT 


due to wind is computed as follows: 


X-COMP = -cos(WDIR-TH) * WKTS * TOF 1.6867 
where: 

X-COMP = down range component of DRT due to wind 

WDIR = true wind direction 

TH = true naedtinns 

WKTS = wind speed in knots 

TOF = time of fall in seconds 

1.6867 = conversion from knots to ft/sec 
This distance added to the PROJECTED DRT gives the WIND 
CORRECTED DRT which is the total DRT the projectile will 


travel in a moving air mass. 


WIND CORRECTED DRT = X-COMP + PROJECTED DRT 


38 





Y 





FIGURE 4. 
NO WIND SOLUTION TO DOWN RANGE TRAVEL 





FIGURE 5. WIND CORRECTED SOLUTION TO DOWN RANGE TRAVEL 
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3. Observed Slant Range 


The experimental data taken from the A-6E aircraft 
during actual said conditions provides many useful delivery 
parameters, sensor readings, and intermediate calculations. 
However, the one piece of information needed to make the 
desired comparison, down range travel, was missing. But 
various other parameters were available, and DRT could be 
reconstructed by several different methods. The problem 
was to determine which method is the most accurate. 

a. Method for DRT Calculation 

Slant range to the target and search radar 
depression angle (look down angle from flight path vector) 
is made available by the search radar. Down range travel 


simply becomes 


DRT = cos(DEPANG) * SR 
where: 
DEPANG = search radar depression angle 


SR = search radar slant range to target 


A second method takes advantage of the aircraft's 
present position altitude and target altitude to compute 


vertical separation and down range travel. 


VERT SEP = PPA - TGT ALT 


DRT = cos arsin(VERT SEP / SR) * SR 
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where: 
VERT SEP = vertical separation 
PPA = inertial derived present position altitude 


TGT ALT = target altitude 


A third and final method of computing down 
range travel uses the vertical separation generated by the 


ballistics program on board the A-6E aircraft. 
DRT = cos arcsin(VERT SEP / SR) * SR 


b. Error Analysis 

An error analysis was conducted to determine 
which method would produce the most accurate value for 
DRT. This would yield a maximum error bound which can be 
expected in DRT due to this method of calculation. The 
third relationship proved to have the smallest error bound 
and was eventually used to reconstruct the down range 
travel from the experimental data. 

For example, the maximum error bound on the 


relationship 
DRT = cos arsin(VERT SEP / SR) * SR 


is the sum of the partial derivatives of DRT multiplied by 


their tolerances. 
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DRT 
VERT SEP 


Apee| = 


Avert SEP| - Boar [Ase 








The third method of computing down range travel will serve 


as an example for this procedure. 
DRT = cos arsin (VERT SEP / SR) * SR 


The partial derivative of DRT with respect to VERT SEP 


becomes: 
© DRT 7 SR 


© VERT SEP (SR)2 - (ALT) 


The partial derivative of DRT with respect to ALT becomes: 


OvRT _ ~ALT 


oO ALT “ ‘sail : (ALT) ~ 


the error bounds are: 


VERT SEP = +0.5 feet (rounded to the nearest foot) 
SR = +5.0 feet (rounded to the nearest 10 feet) 
Therefore, the maximum error bound that can be expected from 


rounding error of actual delivery data becomes: 


Aprt = 0.5(SR) + 5.0(ALT) 


> / CSR) ~ (ALT) 


This analysis was performed using all three 
methods of constructing down range travel. The third 


method was found to be the most accurate with a maximum 
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error bound of 5.2 feet and was the method used to construct 
the down range travel from the freeze data. 
4. Hit Distances 

The hit coordinates of the experimental data were 
not utilized because the serpaeer bases its calculations 
on the location of the search radar cursors. If the cursors 
are not properly placed, the weapon, most likely, will not 
hit the target. However, for the slant range measured by 
the search radar set, the down range travel and time of 
fall will be calculated accurately and the weapon will hit 


the ground in the proximity of the cursor placement. 
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VI. PRESENTATION OF RESULTS 


A. FORTRAN VS. PIM 

Appendix A contains the results of comparing identical 
ballistics algorithms: a FORTRAN program which is a Naval 
Postgraduate School modification of the Naval Weapons Center 
BOEING algorithm, and a PIM version of the same algorithm 
(Ref. 3). The difference between the procedures is that 
the FORTRAN program uses the standard IBM 21 to 24 bit 
mantissa for its floating point number, whereas the PIM 
version uses a 16 bit mantissa. 

A summary of the results is presented in Table 2. 
With the exception of weapon number five (MK-76 MOD-5) 
the largest average difference in down range travel (DRT) 
was 1.6 feet and the maximum absolute difference in DRT 
was 1/ feet, which occurs when DRT is 8,332 feet. Weapon 
number five is suspected to have a coefficient error, 
although none have been discovered. 

Overall these results indicate that the sixteen bit 
mantissa is sufficiently accurate to perform the ballistics 


algorithm. 
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TABLE 2. SUMMARY OF FORTRAN VS. PLM RESULTS 
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B. BALLISTICS TABLES VS. FORTRAN RESULTS 

| The ballistics tables were compared against the results 
of the FORTRAN version of the ballistics algorithm. Appen- 
dix B contains the results for a variety of initial condi- 
tions and weapon types. The results indicate substantial 
discrepancies between the ballistics tables and the FORTRAN 
program. The latest version of the Naval Weapons Center's 
ballistics algorithm (including revised drag coefficients) 
still does not resolve these differences. — Table 3 gives a 
summary of the data in Appendix B. Weapon number five 
(MK-76 MOD-5), which was used in the me cvteeneas data 
shows reasonable accuracy for the range of parameter 


values used in the experiment. 
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AVERAGE - AVERAGE | MAXIMUM 
PER CENT PER CENT PER CENT 


IN DRT IN TOP IN TOF 
(feet) (seconds) (seconds) 


, 
: 
Ts [ease foams [en | ae 


* used invalid drag and mach coefficients 





TABLE 3. SUMMARY OF BALLISTICS TABLES VS. FORTRAN RESULTS 
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C. EXPERIMENTAL DATA VS. FORTRAN ALGORITHM 

Experimental data gathered by bombardiers from Naval 
Air Station Whidbey Island at Boardman bombing range is 
given in Appendix C. Because the ballistics table and 
the FORTRAN algorithm agree reasonably well, the FORTRAN 
algorithm was used in place of the ballistics tables for 
convenience. The down range travel and time of fall were 
calculated from: 

1. The experimental data. 

2. The FORTRAN algorithm with the old drag coefficients. 

3. The FORTRAN algorithm with the new drag coefficients. 
The calculations are described in section (5.B.2, Wind 
Effect), and the results are summarized in Table 4. 

The substantial discrepancy between the experimental 
data and the results of the ballistics algorithm cannot be 
dismissed. Because the ballistics algorithm's down range 
travel agrees with the official Navy ballistics tables to 
within 0.2%, the 12% discrepancy in DRT leads to the con- 
clusion that either the instrumentation on many different 
aircraft indicated erroneous readings, or the behavior of 
the weapons is substantially different from the behavior 
described by the ballistics tables. Additional data under 
more precise initial conditions would have to be gathered 


before any definite conclusions can be drawn. 
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TABLE 4. SUMMARY OF BALLISTICS TABLES VS. ACTUAL 
DELIVERY DATA 





AVERAGE ABSOLUTE 
DIFFERENCE IN 
DOWN RANGE TRAVEL 
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DIFFERENCE IN 
TIME OF FALL 











AVERAGE PER CENT 
DIFFERENCE IN 
DOWN RANGE TRAVEL 







AVERAGE PER CENT 
DIFFERENCE IN 
TIME OF FALL 
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VII. CONCLUSION 


The radical cost reductions in computer hardware 
brought about by large scale integration (LSI) have intro- 
duced an opportunity to construct micro-computer based 
airborne tactical systems which reduce the hardware costs 
by at least an order of magnitude. To establish the 
feasibility of constructing such a system requires that | 
two questions be answered in the affirmative. 

l. Is the computation sufficiently accurate? 

2. Is the computation fast enough to satisfy real 

time requirements? 

This study concentrated on the first question. As the 
results indicate, the 16 bit floating point mantissa is 
sufficiently accurate for the ballistic calculations. 

As a byproduct, the BOEING-Naval Weapons Center algo- 
rithm was compared with the published ballistics tables. 
Although some of the weapons displayed close agreement, 
others revealed substantial discrepancies which remain 
unresolved. 

The most significant and unexpected finding is related 
to the experimental data generated by the bombardiers based 


at Naval Air Station Whidbey Island. If the initial 
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conditions recorded by the aircraft's instruments are used 
to predict where the weapon would impact the ground, then 
Bie ballistics tables predict that the weapons land con- 
sistently more than 10% short of where they actually landed. 
Either the recorded initial conditions are incorrect or the 
ballistics tables do not predict reality for this weapon. 
Although the micro-computer is substantially slower in 
executing arithmetically complex tasks when compared to a 
mini-computer, several micro-computers can be used as 
dedicated machines for specific tasks. Such a distributed 


system can operate sufficiently fast to solve the real time 


problem. 
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APPENDIX A 


This appendix compares the results of the FORTRAN and 
PIM versions of the ballistics algorithm. The absolute 
difference in down range travel and time of fall is 


presented. 
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APPENDIX B 


This appendix compares the NAVAIR O1-1C-1T-1 Ballistics 
Tables with the results of the FORTRAN version of the 
ballistics algorithm. The difference in down range travel 


and time of fall is presented. 
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APPENDIX C 


This appendix is a listing of the output from the cockpit 
of fomious A-65 aircraft recorded at the EER ne ceapon 
was released. 

The following is a brief description of the parameter 


headings as they appear on the listing. 


TAS = true airspeed in knots 
TH = true heading in degrees relative to true north 


WDIR = wind direction in degrees true 


A 


wind speed in knots 

GT = ground track in degrees true 

GS = ground speed in knots 

RA = release angle in degrees 

VZ = vertical velocity in knots x 10 

VSEP = vertical separation in feet 

DA = depression angle (search radar) in degrees 


SR 


slant range to target in feet 


TOF time of fall in seconds 
HIT = his distance from target in feet 


AZ = hit azimuth in clock code 
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APPENDIX D 


Thi aie ean contains the output which os the 
experimental data with the FORTRAN version of the ballistics 
algorithm using the old set of mach, drag, and weapon coef- 
ficients. The weapon used for this experiment was the MK-76 
MOD-5 25 pound practice bomb. 

The coefficients are all assigned in the DECODE sub- 
routine (see Ref. 3 for further explanation of these 
coefficients). Both sets of coefficients, old and new, 


are summarized here for the reader's convenience. 


Weapon Coefficients 


Old Coefficients New Coefficients 
IREF = 2 IREF = 2 

IBOTH = 1 IBOTH = 1 

ITYPE = -l ITYPE = -1l 

DMAX = 3.0 DMAX = 6.0 
CFORM1 = 0.0039077 CFORM1 = .1064453 
CFORM2 = 0.0 CFORM2 = 0.0 
DKG1 = 0.0063648 DKG1 = -.0043918 
DKG2 = 0.0 DKG2 = 0.0 

DM1 = 0.0 DM1 = -.270 

DM2 = 0.0 DM2 = -.270 

VMUZ = 0.0 VMUZ = 0.0 

FN = 0.0 FN = 0.0 

VE = 0.0 VE = 0.0 

SL = 0.0 SL = 0.0 

DS = 0.0 DS = 0.0 
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APPENDIX E 


This appendix compares the experimental data with the 
ballistics algorithm using the new coefficients for drag, 


mach, and weapon type. 
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